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Abstract 

The problem of finding a triangulation of a convex three-dimensional 
polytope with few tetrahedra is ./VP-hard. We discuss other related com- 
plexity results. 

1 Introduction 

A triangulation of a d-dimensional convex polytope P is a set of d-simplices 
whose union is the polytope, their vertices are extreme points of P, and any 
two simplices in it intersect in a common (possibly empty) face. The size of a 
triangulation is the number of its full-dimensional simplices. In this paper we 
discuss the computational complexity of finding small size triangulations of a 
convex polytope. We discuss in particular the case of minimal triangulations, 
i.e. those with smallest possible size. 

This geometric minimization problem arises, for example, in Algebra and 
Mathematical Programming. For example, minimal triangulations of the d-cube 
have been extensively studied (see references in Jll], due to connections 

with the simplicial approximation of fixed points of continuous maps (sec [27j ) . 
Optimal size triangulations appear also in the polyhedral techniques in Alge- 
braic Geometry [ p6| . Understanding minimal triangulations of convex polytopes 
is related to the problem of characterizing the /-vectors of triangulations of balls 
and polytopes (see open problems in Q). In fact, the study of minimal trian- 
gulations of topological balls also received attention due to its connections to 
data structures, in the calculation of rotation distance of binary trees fl25fl . 
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The computational geometry literature has several papers interested in find- 
ing triangulations of optimal size @, |l3). In 1992 Bern and Eppstein asked 
whether there is a polynomial time algorithm to compute a minimal triangula- 
tion of a 3-dimensional convex polytope (open problem 12 in section 3.2 Qj). 
Our main result shows that, under the hypothesis P ^ NP, such an algorithm 
cannot exist: 

Theorem 1.1 Given a convex 3-polytope P with n vertices and a positive in- 
teger K , deciding whether P has a triangulation of size K or less is an NP- 
complete problem. 

We have the following corollaries (the second result was recently obtained in 
|pOf via a direct transformation to 3-SAT): 

Corollary 1.2 1. Finding a minimal- size triangulation of a convex polytope 
of any fixed dimension d > 3 is NP-hard. Clearly, the same holds when 
the dimension is not fixed. 

2. Finding a minimal-size triangulation of the boundary of a convex polytope 
of any fixed dimension d > 4 is NP-hard. Clearly, the same holds when 
the dimension is not fixed. 

Now we discuss the general structure and main ideas of the proof of Theorem 1.1. 

We give a transformation to the Satisfiability (SAT) problem (cf. given 
an instance S of C logical clauses in V boolean variables, is there a truth 
assignment to the variables such that all clauses are simultaneously satisfied? 
We will give a number K and construct a convex 3-polytope, of size polynomial 
in C and V (polynomial size pertains to the binary encoding length), which 
has a triangulation of size at most K if and only if there is a satisfying truth 
assignment. In fact we can restrict our discussion to the special case of the SAT 
problem where each variable appears in three clauses, two of the times negated 
(see page 259 Q). 

Two elementary properties of triangulations will be useful to reach our goal: 
(1) Every boundary triangular facet F of a polytope is contained in exactly 
one tetrahedron of a triangulation. The fourth vertex of that tetrahedron is 
said to triangulate F . (2) Simplices of a triangulation cannot intersect in their 
relative interiors. We will primarily see this behavior in triangles being pierced 
by an edge of the triangulation, a so-called bad intersection. Our proof combines 
techniques presented in the articles l| and p2| : 

Consider a long vertex-edge chain on the boundary of a polytope such that 
the adjacent triangular faces all meet in two points a and b (see Figure Q). The 
proof of the following lemma can be found in Section 2 of jq] : 

Lemma 1.3 Let P be a convex 3-polytope such that the triangles (a, qi, Qi+i) 
and (6, qi, qi+\) for i = 0, . . . , m are among its facets, with the additional restric- 
tion that conv{a, b} n conv{qo, . . . , q m +i} = 0- Let n be the number of vertices 
of P. 
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Then, for each triangulation of P that does not use the (interior) edge (a, b) 
the number of tetrahedra is at least n + m — 3. 




For us, the number n + m — 3 will be a relatively large number, such that 
a triangulation not using (a, b) cannot be small. When using the edge (a, b) 
we can triangulate P c hain = conv{a, b, qo, . . . , q m +i} using the m + 1 tetrahedra 
(a,b,qi,qi+i) for i = 0, . . . ,m. Call Q the (non-convex) polytope we get after 
cutting all these tetrahedra out of P. Let tiq denote the number of vertices of Q. 
Note that n = uq + m. Suppose the number of tetrahedra in any triangulation 
of Q (if there is one at all) is bounded above by some number t. Then we can 
bound the size of a minimal triangulation of P using (a, b) by t + m. Note that 
t does not depend on the length m of the vertex-edge chain. Hence, by choosing 
m large enough (leaving Q as it is) makes 

t + m < uq + 2m — 3, 

and any close-to-minimal triangulation will have to use (a,b). This argument 
still holds when we have many vertex-edge chains of the same length m present 
in other parts of the boundary of the polytope P. If m is large enough, a small 
triangulation is forced to use the edges (a, b) of all these vertex-edge chains. 

We also use the famous non-convex Schdnhardt polytope [0, |l7], [l8|, [ii], H||, |2j| . 
Roughly speaking, a Schonhardt polytope can be obtained by "twisting" the 
top face of a triangular prism in a clockwise direction (see Figure ||) . The three 
quadrangular sides are then broken up and "bent in" , thus creating the non- 
convex (reflex) edges (Bi,Ai + \) that we call diagonals. The resulting polytope 
is non-convex and wc distinguish two triangular faces the bottom (Ai, A2, A$), 
and the top or skylight (Bi, B>2, B3) not having a reflex edge (Note: whenever 
dealing with vertices of a Schonhardt polytope, abusing the notation, by an 
index i + lwe mean (i mod 3) + 1. For example 3 + 1 gives 1). 

The reader can easily verify that Schonhardt polytopes cannot be triangu- 
lated using only its six vertices. Imagine the Schonhardt polytope is glued along 
its bottom face to a bigger polytope. Again, the resulting non-convex polytope 
can only be triangulated if its top triangle (the skylight) is visible by another 
vertex (we will rigorously define visible later, but it does correspond to the in- 
tuition that every point of the skylight can be joined to the vertex by a straight 
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Figure 2: A Schonhardt polytope 



line segment. We will show that the triangular cone defined by the hyperplancs 
containing the faces (£?£, -Bj+i, -Aj+i) contains exactly the points that can view 
the skylight. For this reason we will call this cone the visibility cone. 

Now we convexify the Schonhardt polytope by attaching three circular vertex- 
edge chains opposite to the concavities. This will give us a convex polytope that 
we call the cupola, see Figure ^. The cupola is usually glued along its bottom 
face to a bigger convex polytope and obtain a convex polytope P. We can com- 
bine what we know about vertex-edge chains and about Schonhardt polytopcs. 
Namely, in order to have a small triangulation of P, the three diagonals of the 
Schonhard polytope inside the cupola have to be used. But then, the vertex 




Figure 3: The visibility cone and an invisible vertex 
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triangulating the skylight of the cupola must not be obstructed from seeing the 
skylight by the diagonals. Hence the vertex has to lie in the visibility cone of 
the cupola. 

In [^] Ruppcrt and Scidel used SAT to prove that it is NP-complete to decide 
whether a non-convex polyhedron admits a triangulation. Their constructions 
used Schonhardt polytopes, and in particular their visibility cones, to do the 
transformation. In our case, because we need convexity, we glue cupolas, instead 
of Schonhardt polytopes. They are glued to a bigger frame polytope along their 
bottom faces. Similar to [2^| , we have variable cupolas and clause cupolas. The 
visibility cones of the variable cupolas contain only two truth-setting vertices, 
one for false and one for true. The visibility cones of the clause cupolas contain 
as many literal vertices as there are literals in the logical clause. Each variable 
must choose between a "true" or "false" value. Inside each clause at least 
one variable will be chosen to be true (to satisfy the clause). We model these 
logical choices by the geometric choices of which vertex in the visibility cone of 
a (variable/clause) cupola is used to triangulate the skylight. In addition, our 
polytope satisfies some blocking conditions: the tetrahedron spanned by the top 
face of clause cupola and a literal vertex coming from a negated variable will 
improperly intersect the tetrahedron spanned by the top face of the cupola of 
variable Xi and the truth-setting vertex corresponding to true. In this way the 
choices made for the truth values of the variables and for the literals satisfying 
the clauses will be consistent. We will call our polytope the logical polytope 
because it comes from a logical formula. 

For the polynomial transformation (from SAT) we need to give an algorithm 
to compute the coordinates of the logical polytope. The binary encoding length 
of the polytope, as well as the runtime of the algorithm, have to be polynomial 
in the encoding length of the SAT instance. Each step of the construction will be 
polynomial, this is a delicate point in the formalism of our argument. We apply 
a sequence of these constructions (polynomially many). The coordinates of 
the vertices of the polytope are potentially singly-exponential, but their binary 
encoding length is guaranteed to be polynomial. 
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Elementary steps of construction include operations such as taking the join of 
two or three points, intersecting planes and lines, putting points on polynomial 
curves, etc. The coordinates of the resulting construction elements are therefore 
polynomials in coordinates of the input elements. On the other hand, we will 
have requirements on the positions of the points with respect to some planes 
or other points on lines etc. All these conditions can be formulated as strict 
polynomial inequalities in coordinates of the construction elements. An essential 
element of our construction is that our systems of strict polynomial inequalities 
will depend on one single parameter e. All these polynomial inequalities are 
satisfied at e = 0, but an additional requirement for us is e > 0. The following 
lemma describes a polynomial algorithm to find a number eo such that all < 
e < eo solve the inequality system. 

Lemma 1.4 1. Suppose p(e) = a^e rf + • • • + a\e + ao is a polynomial with 

p(0) > 0. Let e (p) := min [ 1, — - — — — — — — p- ] . Then for < e < 

V 2(|ai|H \-\a>d\)J 

e oip) we have p(e) > 0. 

Hence, the construction of eo can be done in time polynomial in the encod- 
ing length of the coefficients of p, and eo has polynomial encoding length. 

2. pi,...,pi are univariate polynomials such that pi(0) > 0, ...,pi(0) > 
then there is a rational number, eo > ; such that pi(e) > 0, . . . ,pi(e) > 
for all < e < eo . Moreover, the encoding length of eo is bounded by a 
polynomial in the encoding length of the coefficients of ' P\, ■ ■ ■ ,Pi- 

Proof: For < e < 1 we have that a^e' > — |a,|e. The reason is that for eij > 0, 
ait 1 > > — |fflj|e, and for a.j < 0, a^ 1 > a^e = — \di\e. Hence for < e < eo(p) 

d d 

p(e) > ~l a »l e + a o > - N ^/ + a > 0. 

i=l i=l ^ Z^i=l 

For the second part, take the value eo(pi, ■ ■ ■ ,Pr) '■= rnin(eo(pi), . . . , eo(Pr))- 
Now all the conditions are simultaneously satisfied. □ 

Of course, in general the real solutions of a multivariate system of inequalities 
coming imposed by geometric requirements may be empty, but our steps of 
construction reduce everything to sequentially solving easy univariate systems 
of inequalities. 

Here is the organization of our paper. In Section |^ we discuss useful proper- 
ties of Schonhardt polytopes and of cupolas. Later in the section we explain how 
to construct and glue cupolas that have a prescribed visibility cone and how to 
construct visibility cones that fit our purposes. The polynomial transformation 
to SAT is presented in Section ^|. From a given SAT instance we construct a 
frame polytope to which we then glue the cupolas. In the final Section 0, we 
discuss the consequences of our result and related complexity problems. As 
a complement of our main theorem we present a family of polytopes (the so- 



called stacked polytopes) for which the decision of Theorem 1.1 can be solved 
in polynomial time. 
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2 Basic Building Blocks 

We recall the notion of beyond a face (see p8fl): A point p is beyond a face F of a 
polytope P if it (strictly) violates all inequalities defining facets of P containing 
F, but it strictly satisfies all other inequalities that define other facets of P. The 
polytope Pbeyond f is the (closure of the ) set of all points beyond F. We denote 
by P \ F the polyhedron defined by all facet-defining inequalities that do not 
hold with equality in F. This is exactly P U Pbeyond f- In our constructions we 
will often put one or more points beyond some face, and then take the convex 
hull. This will only destroy the facets containing this face, and introduce new 
ones containing the new points. We will say we attach one polytope P to another 
Q along a facets Fp of P and Fq of Q if P C Qbcyond f q an d Q C Pbeyond f p ■ 
It is important to observe that the convex hull of their union contains both the 
face lattices of P and Q without, of course, Fp and Fq. 



2.1 The Schonhardt Polytope 

Let us turn to a well-known example of a non-convex non-triangulablc polytope, 
the so-called Schonhardt polytope (named after its first occurrence in [^3|. Sec 
also fill). For the notion of non-convex polytope and what it means to trian- 
gulate them we refer to |10[ . 

Definition 2.1 A Schonhardt polytope (Figure^) is a non- convex polytope with 
six vertices A\, A 2 , A3, B±, B2, and B3 and facets {A\, A2, A3), {B\, B2, B3), 
{A 1 ,B 1 ,A 2 ), (B 1 ,A 2 ,B 2 ), {A 2 ,B 2l A 3 ), {B 2 ,A 3 ,B 3 ), (A 3 , £3,-81), and 
{B3, B\, A\). At exactly the edges (Bi,A 2 ), (B 2 ,A3), (i?3,Ai) the correspond- 
ing facets are to span an interior angle greater than ir ( the edges are said to be 
reflex,). These edges are called the diagonals of the Schonhardt polytope. The 
top face (.Bi, B 2l B3) is called the skylight of the Schonhardt polytope. 

Six points are said to be in Schonhardt position if they are the vertices of a 
Schonhardt polytope. We say that the skylight is visible from a point x ( or x is 
able to see the skylight, or x is a viewpoint of the skylight) if the tetrahedron 
spanned by x and the skylight does not intersect any of the diagonals in their 
relative interior. The visibility cone of the Schonhardt polytope is the triangular 
cone bounded by the planes B\B 2 A 2 , B 2 B3A 3 , and B3B1A1. See Figure 



r 



The use of the word "skylight" is motivated by the idea that the skylight 
triangle is a glass window and light comes through it illuminating the interior 
of the Schonhardt polytope defining a cone of light. It is obvious that this 
non-convex polytope cannot be triangulated (without adding new points): The 
fourth point of the tetrahedron containing the skylight must be one of A\, A 2 , 
or A3, but the diagonals "obstruct the view" of the skylight from these vertices. 

It is our intention to patch the sides of the Schonhardt polytope with vertex- 
edge chains in order to convexify it (and then glue it to a frame polytope). 



According to Lemma 1.3, a small triangulation of this convex polytope must 
necessarily contain the diagonals. In this case, the fourth point of the tetrahe- 
dron containing the skylight also has to be able to see the skylight. We will show 
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where to place the vertex-edge chains in order for them not to be visible from 
the skylight. Hence, the triangulating vertex has to lie beyond (A\, A2, A3), we 
will show that it has to lie in the visibility cone. 

Lemma 2.2 Let A±, A2, A3, B\, B2, B3 be six points in Schonhardt position. 
We denote by Ca,b the convex hull of the six points. Then 

1. All orientations of simplices spanned by four of these six points are de- 
termined up to one global sign change. As a consequence, the six points 
are in convex position, and their convex hull Ca.b is an octahedron that 
has (Ai, A2, A3) and {B\, B2, B3) as facets and it has edges (Ai,Bi + i) 
(1 = 1,2,3;. 

2. There are no points that can see the skylight (B\, B2, B3) and, at the same 
time, (1) are beyond either of the edges (A, ,_Bj + i) of Ca.b, and (2) are on 
the side of the plane B1A2B3 opposite to B2 or similarly for the analogous 
planes B1A3B2, B2A1B3 and the points B3, B\ respectively. 

3. The visible points beyond the facet {A\, A2, A3) of Ca,b are exactly the 
points that are also in the visibility cone of the Schonhardt polytope. 

In what follows we will use the language of oriented matroids. For the theory 
of oriented matroids we refer to || and |2g] . Here we only sketch the necessary 
definitions and how they are related to the notion of visibility. The orientation 
of a simplex (xi,X2,X3,Xi), is defined as 

r i , . ( Xl X 2 X3 Xi \ 

[x\,X2,X3,x^\ = sign det I 1 ^ 1 1 )• 

All such orientations make up the chirotope of an oriented matroid (see page 
123 in @). 

Given the oriented matroid of points x\, . . . ,x n in d-space, its circuits are 
functions C : {x\, ■ ■ ■ ,x n } 1— ► {+, —,0} that correspond to so-called minimal 
Radon partitions. This means that the convex hulls of C + = {xi\C(xi) = +} 
and C~ = {xi\C(xi) = — } intersect in their relative interiors, and C + and C~ 
are minimal at that. It is easy to check that the function 

omit Xi 



C(x) = \ ' ■ ' ■ i^+i] > iix e {xi,...,x d+1 }, 

otherwise, 

defines a circuit if it is not identical 0. In fact, all circuits can be obtained this 
way. We will compute circuits to use a argument of the following form: x does 
not see the skylight if and only if there is a circuit such that the positive part 
is one of the diagonals and negative part is the set containing x and a subset 
of vertices of the skylight. Since then the tetrahedron spanned by x and the 
skylight is pierced by the diagonal. 
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Important tools to compute simplex orientations are the Grassmann-Pliicker 
relations (see Section 2.4 in For points a, b, xx,...,x& they state that the 
set of signs 



{[a,b,xi,x 2 ] ■ [a,b,x 3 ,Xi],-[a,b,xx,X3\ ■ [a,b,x 2 ,Xi], [a, b, x x , x±] ■ [a,b,x 2 ,x 3 ]} 

is either identical or contains both a + and a — . The typical use of the 
Grassmann-Pliicker relations is to deduce one orientation when the others are 
known. We can read the orientations of some of the different tetrahedra from 
two-dimensional projections (drawings) of the point configurations as in Fig- 
ure ||. We use a left-handed rule system, i.e. we decide whether the triangle 
(xi, x 2 , x 3 ) is oriented counterclockwise (+) or not (— ), also if 2; 4 is on our side 
of the plane spanned by x\, x 2 , and X3 (+) or not (— ), and multiply these two 
signs to obtain the orientation [x\, x 2 , X3, £4]. 



Proof of Lemma 2 A 



1. In a Schonhardt polytope, the simpliccs {A\,A 2 , A3, B{), and (A\, A 2 , A3, B 2 ) 
have the same orientation since edges (Ai,A 2 ) and (A 2 , A3) are both incident 
to facet (A±, A2,A$) and they are both non-reflex edges. 

By the above argument, going around the boundary of a Schonhardt poly- 
tope, keeping in mind which edges are reflex, we can determine the orientation 
of 12 simplices up to one global sign change (there are 12 edges). But there are 
( 4 ) = 15 simplices formed by the vertices of the Schonhardt polytope. The re- 
maining three simplices are (Ai, A 2 , B 2 , B 3 ), (A 2 , A3, B\, B3), {A\, A3, Bi, B 2 ). 
The signs are determined by the following Grassmann-Pliicker relations: For 
(Ai,A 2 ,B 2 ,B 3 ) take a = A 1} b = A 2 , x\ = A3, x 2 = B±, x 3 = B 2 , x 4 = B 3 
(the other two by circular index shift). Then: 

{-?, +,-•-} 2 {+,-} 

the equation forces [A\, A 2 , B 2 , B3] ~ +. From the chirotope information it is 
easy to check that all vertices are in convex position (see description of how to 
read the facets of the convex hull from the chirotope in Chapter 3 of ||), and 
that their convex hull Ca,b is indeed an octahedron. 

2. We will show that if a point x lies beyond ^i-E^ of Ca,b, on the side of 
BxA 2 B 3 opposite to B 2 , then (Bi,A 2 ) and the triangle (B 2 ,B 3 ,x) form a 
minimal Radon partition in the set of vertices A\, A 2 , A3, B\, B 2 , B3, and x, 
hence have an interior point in common. This means x cannot see the skylight. 
For this, we compute the following orientations: 



[B l ,B 2 ,B3,x] = + 

-[A2,B 2 ,B 3 ,x] = + 

~[A 2 ,B l ,B3,x] = - 

-[A 2 ,B 1 ,B 2 ,x] = - 

-[A^B^B^Bs] = - 



since (Bi, B 2 , B 3 ) is a facet of Ca,b \ {A\,B 2 ), 
since (A 2 ,B 2 ,B 3 ) is a facet of Ca,b \ {Ai,B 2 ), 
from the assumption on x, 
from the Grassmann-Piicker relation below, 
from Part (1). 
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The necessary Grassmann-Pliicker relation is the one with a = B\, b = B 2i 
Xi = B 3 , x-i = A\, X3 = A2, and £4 = x such that 

{-•?, 

forces [Bi ,B2,A2,x] = — . 

3. If x is in the visibility cone V, then it is, by part (2) of this lemma, on the 
same side as B 3 with respect to the plane B1A2B2. Hence A2 is on opposite 
side of B 3 with respect to the plane B1B2X. Therefore, the relative interior of 
the convex hull of B\ and A2 lies strictly on one side of the plane B1B2X, and 
the tetrahedron (Bi, B2, B 3 , x) on the other side of this plane. Therefore those 
two point sets cannot have points in common. By symmetry it follows that the 
other two diagonals do not obstruct any point of V from seeing the skylight 
either. 

Assume now that a point x is beyond face (Ai, A2, A 3 ), but outside of V, 
i.e. for instance on the A\ side of the plane B1B2A2. We claim that the pair 
{Bi, A2}, {B2, B 3 , a;} forms a circuit in the oriented matroid of the point config- 
uration of the vertices of Ca,b and x. This means that the triangle (B2, B 3 , x) 
is pierced by the diagonal (B\, A 2 ) in the relative interior, hence x is not visible. 

-[B2,B 3 ,A 2 ,x] = -, since (B 2 ,B 3 ,A 1 ) is a facet of C a ,b \ (Ax, A 2 ,A 3 ), 

+ [£>i, B3, A2, x] = +, from the Grassmann-Pliicker relations below, 

— \B\,B2,A2,x\ = +, from the assumption on x, 

+ [Bi, B 2 , B 3 , x] = -, since (Bi, B 2 , B 3 ) is a facet of Ca.b \ {A\, A 2 , A3), 

-\B l ,B2,B3,A2] = +, from Part (1). 

In this case, we have to apply the Grassmann-Pliicker relations twice to get 
[Bi , B 3 , A2 , x] = +. First we deduce [A\, A2, B 3 , x] = — from the Grassmann- 
Pliicker relation with a = A\, b = A 2 , x\ = A3, x 2 = x, X3 = B2, X4 = B3: 

{-■- -•+•?, + •+} 2 {+,-}■ 

Now we use this orientation to formulate a = A 2 , b = B 3 , x\ = A\, x 2 = B\, 
x 3 = B 2 , X4 = x: 

{+•--•+•?,-•+} 2 {+,-} 
in order to get the desired [A2 , B 3 , B\ , x] = — . □ 

2.2 The Cupola 

Definition 2.3 A polytope C is called a cupola if it has the following properties: 

1. the vertices of C are A\,A2,A 3 , B\,B2,B 3; and q£ 3 (k = 0, . . . , m + 
1, (i,j) E {(1,2), (2,3), (3, 1)}/. where q 1 ^ = A< and q% +1 = B j are 
identified. 
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Figure 5: A cupola as part of a larger convex polytope P 



2. the vertices A\,A2,A^, i?i,£?2,-B3 are in Schonhardt position, and 
(.Ai, A2, A3) (the bottom facet) and (Bi,B2,Bs) (the skylight) are facets 
ofC, 

3. the other facets are (-B,, q]f , q^ 3 +l ) and (Aj , q]f , q l kJrl ) for k = 0, . . . , m+1, 
(i I i)G{(l,2) J (2,3),(3,l)}, 

4- the vertices q]^ 2 (k = 1, . . . ,m) lie on the side of the plane B1A2B3 oppo- 
site to £?2- Similar conditions must hold for q k ' 3 and q^.' 1 . 



Proposition 2.4 Let C be a cupola which is part of a larger polytope P, 
i.e. Q = P \ C is a convex polytope and Q and C share the common facet 
(Ai, A2, As). Let n be the number of vertices of P. 

If T is a triangulation of P with the property that the fourth point of the 
tetrahedron containing the skylight of C is not in the visibility cone of C , then 
there are at least n + m — 3 tetrahedra in the triangulation. 

Proof: If the vertex triangulating the skylight of C is a vertex of on a vertex-edge 



chain of C, then it does not see the skylight by Definition 2.3 (4) and Lemma 



2.2 (2). If it is in Q instead, then it has to be beyond the face (Ai,^,^) of 



C. Hence by Lemma [2.2| (3) it cannot see the skylight either. Therefore the 



triangulation T does not use one of the diagonals. By Lemma 1.3 the number 



of tetrahedra is at least n + m — 3. □ 



Later, Lemma 1.4 will be used to guarantee that we can place a point beyond 
a certain face. 



2.3 Constructing a Cupola from a Visibility Cone 

In this subsection we will show that cupolas can be attached to any face of a 
frame polytope using intermediate polytopes and that the visibility cone can be 
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prescribed. The following theorem does not have the full strength we need for 
the construction. In Section 3, we will use a slightly stronger version which we 
will present at the end of this section. However, this theorem captures the main 
ideas used to construct a cupola. 

Theorem 2.5 (Cupola Construction from a Given Visibility Cone) Let F be a 
facet of a 3-polytope P, and V be a triangular cone such that F(~) V is a triangle 
in the relative interior of F, and m be a positive integer. Then there is an 
m-cupola C beyond F of P such that P is beyond (A±, A\, A3) of C and such 
that V is the visibility cone of C . Moreover, the input length of C is polynomial 
in the input lengths of P, V and m. 

Before we come to the proof, we will exhibit a necessary condition of the 
visibility cone V of a cupola C and the facet the cupola is being glued upon. It 
will imply that we cannot directly attach a cupola to a face (as in p2j), but we 
have to construct an intermediate polytope first. 




Figure 6: Collinearity condition in the base triangle of a cupola 

Lemma 2.6 Let A\, A%, A3, B\, B2, B3 be vertices in Schdnhardt position. 
Define l\ to be the intersection line of planes B3B1A1 and B1B2A2, lines I2 
and I3 are defined accordingly (Figure^, note that they contain the extreme rays 
ofV). The lines l\, I2, and 1 3 intersect the relative interior of the bottom face 
(A\, A2, A3) of a cupola C. The intersection points D\, D2, and D3 are forced 
to have the following collinearities: A\D\D2, A2D2D3, and A3D3D1 . 

Proof: l\ enters the Schonhardt polytope S in point B±, runs along facet 
(Ai, Bi, B3) until it reaches the edge (Ai,B 3 ) where it goes into the interior 
of S. Then the relative interior of (Ai, A2, A3) contains the point D\. In this 
way, Z?i, D 2 , A 2 arc all on the planes A1A2A3 and B1B2A2. □ 



Proof of Theorem 2.5: We proceed in three steps. The lines l\, I2, h are defined 
as in Lemma o. 
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The bottom triangle (A\, A 2 , A3). We will now construct an intermediate poly- 
tope beyond F which will have a triangular facet (Ai,A 2 , A3) which is (1) par- 
allel to F, and which is (2) intersected by the cone V in a triangle (Di, D 2 , D 3 ) 
in the relative interior such that (3) the collinearity condition from Lemma 2.6 
holds. 

To do this, we place a plane H parallel to and slightly above F such that 
the intersection points Di of H and U (i = 1, 2, 3). Also H has to be so close to 
F that the U do not cross between H and F. By prolonging the line segment 
D3D1 slightly beyond D\ (staying in P \ F) wc obtain point A\, analogously 
construct A 2 and A3 (Figure 0). Taking the convex hull of F and the points A\, 
A 2 and A3 gives then the intermediate polytope, whose face (Ai, A 2 , A3) has 
the collinearity condition. These constructions are polynomially constructible 
in the sense of Lemma 1.4. 




D, 



A 1 ^D 1 D 2 A 2 



F 



V 




Figure 7: Building the intermediate polytope for the cupola 

The frame of the cupola. As in the construction of the bottom facet (A\, A 2 , A3), 
we place a plane H' parallel and slightly above this facet. The intersection of 
H' and V is the triangle (B±, B 2 ,B3) (B\ is on the same extreme ray of V as 
D\ and so on). See figure |[ 

It is clear from the construction that triangles (D\, D 2 , D3), (^1,^2,^.3), 
and (Bi, B 2 ,B3) are parallel and all oriented the same way. Therefore it is not 
hard to check that the points Ai, A 2 , A3, B\, B 2 , and B3 are vertices of a 
Schonhardt polytope whose visibility cone is V. Polynomiality of this part of 



the construction follows from Lemma 1.4 as well 



Attaching the vertex-edge chains. Now that the frame of a cupola done, 
i.e. the vertices A\ , . . . , B3 in Schonhardt position, it remains to patch the 
key structures of Lemma L3, the vertex-edge chains qi' (i — l,...,m, 



(J,k) <E {(1,2), (2,3), (3,1)}), to the sides of the frame conv(Pl) {A u . ..B 3 }). 

Given triangular faces (a, go, Qm+i) and (b, qo, q m +\) of a convex polytope P 
and a plane G which (strictly) separates points qo and q m +i- We claim that 
we can construct points q±, . . . , q m beyond the edge (qo, q m +i) of P such that 



the convex hull of P U {qt} has the properties of Lemma |L3j and such that the 
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Figure 8: Building the frame of a cupola 



points qi, . . . , q m lie on the same side of G as go- Moreover, the input length of 
the constructed points is polynomially bounded in the input length of P and G. 

By applying our claim three times, we will conclude our proof. The vertices 
q?' J+ are placed beyond edge (Aj,Bj + i), vertices Bj and Aj+i take the roles 
of a and 6, G is the plane spanned by Bj, A j+i and Bj + 2- It is easy to check 



that this is exactly what we want for Lemma 1.3 and for the cupola conditions. 

Now we prove the claim. We will put the points qi (i = 1, . . . , m) on a 
parabola segment, beyond the edge (go,g m +i)- Let H be a plane containing go 
and q m +i which also intersects the interior of P. This plane has the property 
that it contains points beyond edge (qo, q m +i)- It is constructible in polynomial 
time. (Let v be the sum of the two normal vectors of planes aqoq m +\ and 
bqaQm+i, and H the plane containing go an d g m +i parallel to v.) 




Figure 9: Construction of the vertex-edge chain 
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Let now D be the intersection point of G and (go, Qm+i)- Let w be a vector 
of direction of the intersection line of G and H , such that starting at D it 
is pointing out of P. Now let E e = D + ew for e > to be specified later. 
For small e, £7 S is beyond (go,g TO +i). Hence the parabola defined according to 
Lemma 2.7, stated and proved below, by p(0) = go, p(l/2) = £/ € , and p(l) = 
q m +i lies entirely in H, and for arguments between and 1 passes just beyond 
(90)9m+i)- Let gj = p{i/(Amj) for i = 1, . ..,m. For small e all those points 
are beyond (go, q m +i) and on the same side of G as go (polynomial conditions, 
use Lemma |l.4[ ). Also, they are in convex position such that the convex hull of 
P U {gi, . . . , g m } has exactly the required face lattice. □ 



Lemma 2.7 Let po, p±, pi three non-collinear points in Mr and to, t±, ti three 
distinct real numbers. Then there is a unique curve p : M — > M 3 such that 
Pa = p(to), pi = p(t±), and p2 = p(^) which is quadratic in every coordinate. 
Furthermore, all points on p(t) are in the plane spanned by po, p\, andp2, and 
they are in convex position. Also a plane containing p{r) and p(l) for some 
r 7^ I which does not contain all of p has all points between I and r on one of 
its sides and all other points on the other side. 



Proof: Since po, pi, P2 have to be on the t$, t±, ti positions of the curve 
p{t)- 



a x + b x t + c x t 2 
a y + b y t + c y t 2 
a, + bd + c,f 2 



we have the condition 



a,. 




1 1 1 

to tl t2 
y-2 f2 4-2 



\ . . \ 

PO Pi P2 



\ 



J 



By the non-singularity of the Vandermonde matrices, there is a unique solution 
to a., & , c. given the pi and t\. 

The curves which arc quadratic in every coordinate are linear transforms 
of the moment curve m(t) = (l,t,t ). This curve lies entirely in the x = 1 
plane, is convex, and has the condition that it intersected by each plane at most 
twice (or it is in this plane). All these properties are invariant under linear 
transformations. □ 



Proposition 2.4 stated that we get a large triangulation if wc triangulate the 
skylight of a cupola by a vertex outside the visibility cone. Now we want to 
estimate how much smaller a triangulation is if we use a vertex v in the visibility 
cone instead. Wc give a relatively small triangulation of the cupola and of the 
space between the bottom face (A\, A2, A3) of the cupola and the triangular 
face F of P with the help of the vertex v. 
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Proposition 2.8 Let F be triangular face of a polytope P, and C an m-cupola 
attached to it according to Lemma 2.5. Let v be a vertex of P in the visibility 
cone of C . Then there is a triangulation of conv({v}, F, C) with at most 3m + 16 
tetrahedra. 



Proof: First of all, wc triangulate along the vertex-edge chains using the tetra- 



, to. 



hedra (Bi, A i+ i, q£ , ) for % = 1, 2, 3, and k = 0, 

After removing these tetrahedra, we are left with the union of the Schonhardt 
polytope on the vertices A\, A2, A 3l B\, B2, -B3, and the convex polytope 
conv({i>}, F, (Ai, A2, A3)). This is a non-convex polytope with all edges, except 



the diagonals, being convex (easy conclusion from Lemma 2.2 and the construc- 
tion). Since the specified vertex v is inside the visibility cone, it sees all facets 
of this polytope, except the three facets it is incident to, from the interior. In 
particular, we can form tetrahedra of all these facets and v and none of them 
intersect badly. They are 7 tetrahedra for the facets of the Schonhardt poly- 
tope (since we do not count the bottom face) and at most 6 for the rest (the 
convex hull of F and (A\, A2, A3) has — by a planar graph argument — at most 
2-6 — 4 = 8 facets, subtracting 2 for F and (Ai, A 2 , A 3 ) gives 6). □ 



It is this 3to in contrast to the 4m in Proposition 2.4 which makes this 
triangulation better for large to. 



2.4 Constructing a Visibility Cone 

In order to use the cupola as a basic building block, we need to have a visibility 
cone that contains a specified set of vertices and intersects the relative interior 
of some face. Once we have that we can construct the cupola as described in the 
previous section. The set will consist of all vertices lying in a specified plane. 

Lemma 2.9 Let H be a plane which intersects the relative interior of some face 
F of a polytope P , and let S = {vi, . . . , v s } be the set of vertices of P lying in 
H , not including the vertices of F . Let S' = {wi, . . . , uv} a set of points in 
relint(F) D H. It is possible to construct a triangular cone V which intersects 
F in a triangle that lies in the relative interior of F and V contains S and S' 
in its interior and no other vertex of P. 

The reader may not see at this point the purpose of the set S' , but we will 
justify it at the end of this section. 

Proof: P fl H is a polygon. Without loss of generality, F n H is horizontal and 
situated on the top of the polygon P n H (see Figure |n|). Let I be the line 
connecting the leftmost point of S' and leftmost vertex of S (the one encoun- 
tered first when walking around P <1 H counterclockwise, starting at F n H). 
Analogously, let r be the line connecting M and the rightmost vertex of S. 

The area between I and r (in H) is already a cone containing S and no other 
vertices of P. We will perturb it in a way that the other conditions are satisfied 
as well. 
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First shift I and r parallely outwards, guaranteeing that they still intersect 
F fl H in its relative interior (easy open conditions); we obtain V and r'. Also, 
let /' be a line in H parallel to F just outside P, i.e. such that V and r' intersect 



/' in the same order as F n H (again using Lemma 1.4). 





Figure 10: Construction of I and r, then I' and r' (viewed in H) 

Now we will rotate H about V and r' and /', getting three planes bounding 
the desired triangular cone: Let H be oriented in some way, and anx > bn be 
its defining inequality. Let v be some point which lies on the positive side of 
H . Let Gi' be the plane through I' and v. By construction, all vertices in S lie 
on the same side of Gy , so we can orient it such that S is on its positive side. 
Let ay x > by be its defining inequality. Perform the same construction for r' 
and /' obtaining G r > and Gfi, also orienting them in a way that v is on their 
respective positive sides. Let G\, be the plane defined by (a H +eay)x > bn+£by. 
This plane contains I' and for small e it is very close to H . Hence, it is the 
rotation of H about I' in the direction of plane Gy . Also let G £ r , be defined by 
(an + ea r i)x > bn + eb r i, and G £ f, be defined by (—an + £df>)x > —bn + £&/'• 

Obviously, all points in S and in S' are on the positive sides of the planes 
Gf, , Gp, and Gp. For small e > 0, these planes do not "sweep" over vertices of 
P which are not in S, and it is easy to see that in this case, there are no vertices 
of P that satisfy all three new inequalities. Also for small e, the points in F 
satisfying all three inequalities define a triangle in the relative interior of F with 
endpoints G\, n G%, n F, G%, n G e fl n F, and G e fl n G e v n F. Hence, the set of all 
points satisfying the three inequalities is a triangular cone V with the desired 
properties. The conditions on e are open polynomial conditions according to 
Lemma 1.4. □ 



This lemma can be used to build one cupola over the facet F . However, there 
might be problems if we keep on constructing around the polytope, like adding 
more cupolas over other facets of P. The visibility cone we just constructed 
might "catch" points we construct later. But these constructions all happen 
beyond facets of P, so we can use the following lemma to construct all cupolas 
one after the other without their visibility cones catching extra vertices. 



Lemma 2.10 Let H±, . . . , H n hyperplanes, intersecting facets F\, . . . ,F n of a 
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F 



Gi, 



G\ 



Figure 11: Rotated hyperplanes, viewed by their intersections with P 

polytope P with the restriction that Fi r\Hj = for all i =/= j. Then Pbeyond F t H 
Hj = for all i ^ j. 

Proof: Assume there is a point u in Pbeyond F t H Hj (i ^ j). Then this point 
also lies in (P \ Fi) (~1 Hj, but on the non-positive side of Pj. Let v be a point 
in Pj n Hj, Then u is also in (P \ P) n Hj (since P,- C P C P \ P t ), but on the 
positive side of Pi. Hence, there must be a point to on the line segment [it, v] 
which is on the hyperplanc containing Pj. The whole segment lies in P \ Pj, 
hence every point on it has to satisfy all of P's defining inequalities except that 
of Fi. So w lies in the facet Pj. But it also lies in Pj (the whole line segment 
does), which contradicts the assumption Pj n Hj = 0. □ 

In Section 3.1 we will need an additional condition: Given a set of lines in 



the plane H (of Lemma 2.G) that pierce the face P, we want to be sure that 
these lines also pierce the skylight of the constructed cupola. (This condition 
will play an important role when we want to force so-called blocking conditions, 
see Section 3: At some point two tetrahedron spanned by two skylights and 
two respective visible vertices v and v' will have to intersect in their interiors. 
This is already guaranteed if the corresponding lines g and g' intersect inside 
the polytope.) 

The next theorem specifies the way in which we will use all the preceding 
lemmas in our construction in Section 3: 

Theorem 2.11 (Full-Strength Cupola Construction) Let H be planes that in- 
tersect facets Fi of a polytope P in their relative interiors such that Fi D Hj = 
for all i + j. Let S { = {v\, ...,<}:= (veH(P)nH i )\F il and L, = {g\, . . . , g\,} 

sets of lines. Assume further that each of the lines g^ pierces the relative interior 
of Fi and is incident to some v\ ■ 

Then we can sequentially construct all cupolas Ci beyond the faces Pj such 
that in the resulting polytope their visibility cones contain Si and no other ver- 
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tices. In addition, the skylight of the cupola over each Fi is pierced by the lines 
in Lj . 



Proof: The theorem follows from the ideas in Lemmas 2.5 and [2.9| . In the 
construction of the visibility cone over facet Fi, we invoke Lemma |2.9| with the 
polytope -PU(J Pbeyond f } ■ The set S[ is of course {lC\Fi\l £ Li}. The cupola 
construction was such that the cupolas over Fj were always beyond the facet Fj , 
so the constructed visibility cones contain no vertices of the other cupolas. In 
order to have the lines in Li pierce the skylight of cupola i we have to alter the 



construction of the cupola in Lemma 2.5: when we put the planes parallel to 



F, we do it in such a way that the triangles (Ai, Ai, A-$) and then (B\, P>2, B3) 
are pierced by these lines. These are both open conditions on the distance of 
the planes to F. □ 



3 The Transformation from SAT 

It is our intention to model the well-known satisfiability problem (SAT) using 
the visibility cones of cupola polytopes. Just as Ruppcrt and Seidel did in ]22] ], 
we will restrict our attention to special SAT instances: each variable appears 
exactly three times, twice unncgatcd and once negated. This is not really neces- 
sary, but simplifies explanations. For our purpose this restriction is appropriate 
because the SAT problem remains NP-complete even for instances where each 
variable or its negation appear at most three times (see references on page 259 
in 10 ) . In addition, note that a change of variables can be used to change a 
non-negated variable into a negated variable if necessary. Also note that if a 
variable appears only negated or only positive the variable and the clauses that 
contain it can be discarded. Finally, if a variable appears exactly once positive 
and exactly once negated then it can be eliminated by combining the two clauses 
that contain the two variables into one. From now on, all logical formulas will 
have the properties that each variable occurs exactly two times unnegated and 
exactly once negated. The formula 

/ = (X 1 V -.X 2 V A 3 V -.Xi) A (-.Xi VI 2 V -A 3 V A 4 ) A (Xi V A 2 V X 3 V X 4 ) 

is such a special SAT formula. The figures in this section will correspond to this 
particular instance. 



In Section 3.1, we will define the logical polytope associated to a given logical 



formula. In Section 3.2 we will compute a number K and sec that the logical 



polytope has a triangulation with < K tetrahedra if and only if the logical for- 



mula is satisfiable. Finally, in Section 3.3 we will give an algorithm to construct 



explicit coordinates of a logical polytope. 
3.1 The Logical Polytope 

In the logical polytope there will be a cupola for each clause and one for each 
variable and its negation. The cupolas will be glued to a frame polytope which 
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resembles a wedge. Look carefully at Figure [12] for an example of the overall 
structure. 




3 clauses 

Figure 12: Sketch of the logical polytope 



Figure |13| displays the view of the lower hull of the frame polytope, seen 
from the inside, i.e. from above. The sharp part of the wedge consists of 2C + 1 
vertices (where C is the number of clauses) cq, . . . , c-ic- We call this part of the 
frame polytope the spine. We attach the clause cupola associated with clause i 
to the triangle (c2j, Cfci+i, £24+2) (shaded in the picture). 





Figure 13: The spine of the wedge: here the clause cupolas are attached 

On top of this wedge structure we will put a series of roofs. They are 
triangular prisms, spanned by the two triangles (x l T , x l F , z l A ) and {x\, z R , z l B ), 
one for every variable Xi of the logical formula. The variable cupolas will be 
attached to the triangular facet (z l L , z R , z B ), the back gables (the triangular faces 
are shaded in Figure |l4| ). 

The variable cupola of variable Xi is such that its visibility cone contains 
exactly the front vertices vertices z % T and z % F . We will use these cupolas to 
read from a small triangulation of the polytope the logical value of variables 
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Figure 14: The roofs, back gables shaded 



according with the following rule: if the truth-setting vertex z l T associated to 
the ith logical variable is used to triangulate the skylight of the cupola for 
variable i, then we set Xi = true. If the truth-setting vertex used to triangulate 
the skylight of the cupola for variable i is instead z l F then Xi = false. 

Beyond the quadrilateral face containing z l T we will place the literal vertices 
x\ and x\ which corresponds to the positive occurrences of Xi in the logical 
formula. Beyond the other quadrilateral face we will place the other literal 
vertex x\ which correspond to the negated occurrence of X{. These vertices are 
in the visibility cones of the three cupolas of the clause where variable Xi or its 
negation appears. 




Figure 15: A roof, back gable shaded, ^-coordinate superelevated 

We list the five conditions on the logical polytopc which are necessary for 
the transformation to work in both ways, i.e. a small triangulation yields a 
satisfying truth assignment for our logical formula and vice versa. 

Definition 3.1 For a given logical formula, a logical polytope is a three- 
dimensional polytope P that satisfies the following conditions: 
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1. (Convexity) The logical polytope must be convex and the face lattice is as 
we just described it. 

2. (Visibility) The literal vertices x\, x\, and x\ are vertices in the visibility 
cone associated to their respective clause cupolas, but of no other clause 
visibility cone. The vertices the only vertices in the visibility 
cones of the ith variable cupola. 

3. (Blocking) This constraint ensures that the assignment of true or false 
values for variables is done consistently, i. e. the positive ( negative ) literals 
can be used to make their clauses true if and only if the variable is set true 
(false). 

Concretely, the tetrahedron spanned by z F and the skylight of the cupola 
of variable Xi intersects the interior of tetrahedron spanned by x\ (by x\) 
and the skylight of the clause cupola corresponding to x\ (to x\). Also the 
tetrahedron spanned by z T and the skylight of the cupola of variable Xi 
intersects the interior of the tetrahedron spanned by x\ and the skylight of 
the clause cupola corresponding to it. See Figure \Tq for an example. 




skylight of clause j 



Figure 16: Blocking for consistent logical values 

4- (Non-blocking) Using the vertex z % T to triangulate the interior of the 
i-th variable cupola should not prevent the non-negated literal vertices 
from seeing their associated cupolas. Concretely, if j is the clause cor- 
responding to the literal vertex x\, then tetrahedra (z l T , z % L , z l R , z l B ) and 
(x\, C2j-2, c 2 j-i, C2j) do not intersect at all. The canonical analogue shall 
hold for x\ and x\ (for x\ replace zt by zf). 
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5. (Sweeping) We follow the same sweeping procedure proposed by Ruppert 
and Seidel jllj/ we will need that 

(a) the variable x\ is to the "left" (negative x direction) of the planes 
C2k-iC2kZp, c 2 kC2k+iZp, and c 2k -iC2k+iz F for < k < C - 1. 

(b) x\ is to the "left" of the planes C2k-iC2k%\, c-2kC2k+ix\, and 
C2k-iC2k+ix\ for < k < C - 1. 

(c) is to the "left" of the planes C2k~iC2kZ F , C2kC2k+\zp, and 
C2k-iC2k+iz F for < k < C - 1. 

(d) z l T is to the "left" j>f the planes C2k~iC2kX 2 , c 2 fcC 2fe+1 a;| ) c 2 k-\C2k+ix\, 
C2k-iC2kx\, c 2k C2k+ix\, and c 2 k-\C2k+\x\ for < k < C - 1. 



3.2 Using the Logical Polytope 

Lemma 3.2 Let P be the Logical polytope, m the number of vertices on each 
vertex-edge chain, n be the total number of its vertices. For a SAT formula 
containing C clauses on V vertices there are polynomials K(C, V) and m(C 7 V) 
such that a logical polytope with m = m(C, V) vertices on each vertex-edge chain 
admits a triangulation with < K = K(C, V) — n + m — 4 tetrahedra if and only 
if there is a satisfying truth assignment to the variables of the logical formula. 

Proof: If a triangulation T of the polytope has < n + m — 4 tetrahedra, then 
by Proposition p.4j the skylight of each cupola is triangulated by a vertex in 
the visibility cone of the cupola. In particular, one of z F and z % T is chosen to 
triangulate the cupola corresponding to variable X; for each i. We claim that 
assigning to Xi the truth value according to this choice (z F ^ false, z F > true) 
satisfies all clauses of the formula. 

Each clause cupola skylight is triangulated by one of the literal vertices, say 
clause j by the positive literal vertex x\ (or x\). By the blocking conditions, it 
cannot be the case that the variable skylight of Xi is triangulated by z F since 
these tetrahedra would intersect badly. So we had set Xi to true. Having x\ 
(or x\) in clause j's visibility cone meant that variable Xi appears unnegated in 
this clause. If the skylight of clause cupola j is triangulated by x\ , by the same 
argument Xi was set to false, and clause j satisfied by the literal -iX;. Hence 
all clauses are satisfied. 

Now we need to prove the converse. If there is a true- false assignment 
that satisfies all logical clauses we must find a triangulation that has no more 
than K tetrahedra. For that we construct a "small" triangulation. There are 
three different kinds of tetrahedra: the ones triangulating the cupolas, the ones 
triangulating the roofs, and the ones triangulating of the rest of the wedge. We 
know how to triangulate a cupola if we know a vertex in its visibility cone (see 



the proof of Prop. 2.£). For the rest we will now follow a sweeping procedure 
which was first described by Ruppert and Seidel p2| . 

The sweeping triangulation proceeds by triangulating "slices" that corre- 
spond to the different variables Xi to Xy, i.e. from right to left. The variable 
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roofs are arranged sequentially for exactly this purpose. A slice is roughly speak- 
ing the part of the tetrahedra between a roof and vertices of the spine. After 
the ith step of the process the partial triangulation will have triangulated the 
first i slices. The part of the boundary of the partial triangulation that is inside 
the logical polytope will form a triangulated disk. We will call it the inter- 
face following the convention of Ruppcrt and Scidcl. It contains the following 
triangles: 



(z T ,c 2 c,zl) and 



(Z T , C2j-2, C2j) 
{zip, C2j-2, c 2j-l) 

and (zp,C2j-i,c 2 j) 



if clause j is satisfied by one of 
the first i variables, or 

otherwise, 



for all j = 1, . . . , C . Before the first step, the partial triangulation is empty. 




ce = c 2 c 



Figure 17: The interface after step 2 



After the last step the partial triangulation will cover the whole logical polytope. 
In general, the vertices of the ith roof will see all triangles of the interface and 
will be used as apexes to form new tetrahedra to add to the current partial 
triangulation. This way the interface will slowly move from right to left. 

Now we describe in detail the triangulation step for the ith variable Xi. 
Since we are only concerned with roof vertices in roof i, we will drop all super- 
scripts. The triangulation step depends on whether Xi is set true or false in 
the satisfying assignment. Let us consider first the case Xi = true: 

The point zt is used to triangulate the interior of the variable cupola as- 
sociated to Xi according to Prop. 2.8. From zt we also form six tetrahedra 



with the following triangles: (z L ,x 3 ,z B ), (x 3 ,z Bl z A ), {z b ,Za,x 2 ), (z B ,x 2 ,Zr), 
(za,xi,x 2 ), and (xi,z A ,z F ). 

Now we come to the part of the triangulation which gave the sweeping pro- 
cedure its name. We form the tetrahedra between x\ and the current interface 
triangles. This is possible by part (a) of condition 5. We also use the tetrahe- 
dron (xi, zt, Co, zp). This is illustrated in the transition from a. to b. in Figure 
|Tj^ . The interface advances to x\, i.e. if (zp, Cj, C&) was an interface triangle 
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Figure 18: Removing the tetrahedra spanned by zt and the shaded triangles 
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Figure 19: The sweep 
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before, now (xi, Cj, Ck) is an interface triangle. Also (zf, C2Cj Z J?) is replaced by 
the triangle (xi, c 2 c, z R ). 

Since is set to true we can use x± to triangulate its clause cupola according 
to Prop. 2.8. We only do this if the clause cupola has not been previously 
triangulated using an other literal vertex. The condition 2 ensures that xi is 
in the visibility cone of the clause cupola coming from the clause that contains 
the unncgated literal Xi. Furthermore, condition 4 ensures that we can actually 
perform this triangulation of the clause cupola without badly intersecting the 
tetrahedra of the variable cupola. In Figure [T^.c. we see that if xi is used to 
triangulc clause j's cupola, then the interface triangle (xi, C2j-2, C2j) is replaced 
by the two triangles (xi, c 2 j-2, czj-i) and (xi, c 2 j-i, c 2 j). 

We repeat this procedure with x 2 , i.e. form tetrahedra with x 2 and the 
current interface triangles, and then use x 2 to triangulate its clause cupola if 
necessary (Figure ^9[d.). This is possible by part (b) of condition 5. We continue 
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by forming tetrahedra using zt as apex (Figure 19 , e, possible by condition 5, 
part (d)). At last, we will include the triangle (c2C, zl,zr)- All these triangles 
are visible by part (d) of the sixth condition. After all these tetrahedra are 
added the interface is ready for the next variable. 




Figure 20: The sweep for Xi = false 



Let us now consider the triangulation step in the case X t is set to be false: 
We use the vertex zp to triangulate the variable cupola as well as seven faces 
of the roof (sec Figure |o|): (z t ,X3,za), (x$,z a ,zb), (x^,z l ,z b ), (z b ,z a ,x 2 ), 
{zb, X2, zr), {za,X2,x\), (x2,x\,zr). The reader can see that on the roof we 
are leaving only the vertex x-$. Next the tetrahedron (zf, z l , zr, C2c) is cut out. 
Hereby the interface triangle (zf, zr,C2c) is replaced by (zf,zl,C2c) (Figure 
^.c). Then IF-j will be used as apex with the triangles in the interface. If 
the negated literal Xi is used to satisfy its clause j, the jth clause cupola is 
triangulated by X3. The interface advances as in the true-case. Then zt can 
be used to form tetrahedra with the triangles in the interface. In the end the 
interface is again ready for the next variable. 

How may tetrahedra can such a triangulation have? Triangulating all cupo- 
las with a vertex in their visibility cones yields at most (3m + 16) (C + V) 
tetrahedra (Prop. |]^). In one step of the sweeping triangulation the tops of 
the roofs are each triangulated using six or seven tetrahedra (if the variable 
is unnegated or negated, resp.). Furthermore, the interface is triangulated by 
some vertices three times (in the positive case by x\, by x l 2 , and by z^) or two 
times (in the negative case by x l 3 and by z % T ). The interface contains in each 
step between C and 2C triangles. Eventually, in cither case there is one more 
tetrahedron (see above). An upper bound for the size of this triangulation is 
therefore 

#T < (3m+16)(C + V) + 7V + 3CV +1 
= m(3C + W) + 16C + 2W + 3CV + 1 

Pt(V,C) 

What is the number of the vertices of the logical polytope in terms of the 
number of clauses and variables? We have V logical variables and C clauses in 
the SAT instance. Say we have m interior points each of the vertex-edge chains 
we added (later we will replace the m by a function of V and C). We observe 
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that we have 3m + 6 vertices in each cupola, hence we have (3m + 6)(V + C) 
for all cupolas. We have in each roof nine vertices, two of them are shared with 
the subsequent roof except for the last roof. Hence the total number of vertices 
in roofs is 7V + 2. We have left only the 2C + 1 vertices along the spine. In 
conclusion, the number of vertices of P is 

n = (3m + 6)(V + C) +7V + 2 + 2C +1 
= m{3C + W) + 8(7 + 137+3 

Pn(C,V) 

We had said before that a "bad" triangulation (where at least one cupola 
skylight is triangulated by a vertex not lying in its visibility cone) has at least 
n + m — 3 = m(3C + 3V + 1) + p n (C, V) tetrahedra. On the other hand a 
"good" triangulation has at most m(3C + 3V) + pt(C, V) tetrahedra. We can 
now set m > p n (C, V) - p T (C, V) and K = m(3C + W) + p T (C, V). Then, if 
a good triangulation exists, its size is smaller than or equal to K, and if not, 
all triangulations are larger than K. Note finally that the chosen m and K are 
polynomial in C and V. □ 



3.3 Constructing the Logical Polytope 

Lemma 3.3 There is a polynomial algorithm that, given a logical formula on 
V logical variables and C logical clauses, produces a convex three-dimensional 



logical polytope as defined in Section 3.1. The coordinates of the vertices of the 



constructed polytope have binary encoding length polynomial in V and C . 

Proof: The construction will be carried out in five stages. By the time we end 
the construction all five requirements of the definition of the logical polytope 
must be satisfied, but three of the conditions will not be met until the last stage. 

1. Give coordinates of the basic wedge, with rectangular faces on top for each 
variable. 

2. Attach the roofs for each variable, giving preliminary coordinates for the 
literal vertices and preliminary coordinates for the points on the lower 
edge (the spine of the wedge) . 

3. Perturb the literal vertices to their final positions. 

4. Perturb the vertices on the spine of the wedge. 

5. Attaching the variable cupolas following the procedures of Section 2. 

In every step we will build a construction element (a point, a line, or a plane) 
whose coordinates are polynomials in the construction elements up to that par- 
ticular moment. Hence, the encoding length of each new construction element 
is bounded by a linear function of the encoding length of the construction so far. 
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The number of construction steps is polynomially bounded in C and V. Hence 
the encoding length of the whole construction is also polynomially bounded in 
C and V. Note however, that the coordinates themselves will in general be 
exponentially large. 

Instead of writing explicit (and highly cumbersome) coordinates for the con- 
struction elements, we rely on Lemma |1.4| to ensure that such coordinates can 
be found if one has really the desire to see a particular logical polytope. A key 
property of Stages 2-4 in the construction is that the geometric conditions we 
want to determine a finite coll ecti on of strict polynomial inequalities in a sin- 
gle variable. Then, by Lemma LA , we know there is an appropiate polynomial 
size solution. In subsequent stages of the construction similar new systems, for 
other independent parameters, will be solved, preserving what we had so far, 
but building up new properties. 

Stage 1: The basic wedge. Consider the triangular prism which is the convex 



hull of the six points c = (0,0,0), C2C = (0, 1,0), 



(0,0,1), 4 = (1,0,1), 



(0, 1, 1), and z R = (1, 1, 1). Sec Figure gT](a). In order to obtain a convex 
structure on the top of the wedge, we consider the function /(x) = x(l — x) + 1. 
The vertices of each roof boundary (that is z l T and z l F as well as z l R and z l L ) will 
face z = f(x). More specifically, z l F = zl^ 1 = (i/V, 0, f(i/V)) and 
(i/V, 1, f{i/V)) for i = 0, . . . ,n. By the concavity of /, the points 
are indeed in convex position and their convex hull, the wedge has the desired 
face lattice (see Figure |2l](b)). 



vt 1 





Figure 21: Construction of the wedge 



So far none of the conditions we want are satisfied (not even partially) . 
Stage 2: The roofs. We will first attach the points z\ and z l B to the quadrilateral 
face (z l L , z l R , z?p, z F ). Then we give preliminary coordinates to the literal vertices 
and to the vertices on the spine. 

Let z\ = l/2-(4+4) + (0,l/3,f roof ) and z B = l/2-(4 + 4) + (0, 2/3, i roof ) 
where t IOO [ is a non- negative parameter that is called the roof height. That is the 
points have the same x coordinate as the midpoint between z l T and z F , y coordi- 
nate 1/3 and 2/3 respectively, and height t roo f over the face (4, z F , z l L ,z R ). Wc 
want to choose t roo { in a way that z\ and z l B are beyond the facet (4, z Fl z l L , z R ) 
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(see Figure |2l](c)). We can easily achieve this by the technique presented in 
Lemma 1.4: The only possibly concave edges are the (zlj,, z\). One restriction is 



therefore that all determinants det(z^,z l L ,z l 7' , z^) have to be positive. These 
arc finitely many open quadratic conditions on t roo f- For i roo f = the points 
z\ an z l B are inside the facets (z\<, z l F , z % L , z l R ), hence the edges in question are 
trivially convex. Wc will get more polynomial constraints on < roo f below and 
then solve all simultaneously to find the suitable roof height. 

The spine of the wedge is still a line. We now put preliminary points 
Co, ... , C2C on this line. Let 

1 3 



U{J) = 2 2C 



and Cj = (0,u(j),0) for j 



0, . . . , 2C - 1, and c 2C = (0, 1, 0) (see Figure |2J). 
As an auxiliary point, let b[ be the barycenter of the points C21-2, c%i-\, and C21 
(1 = 1,..., C). At this moment, this point bi = C2j_i. Later, as we perturb the 
spine vertices 6/ will move accordingly, always bi = l/3(c2z_2 + c%x-i + 



Cl C3 C5 C 7 
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Figure 22: Preliminary coordinates for the spine vertices 



Now we want to give initial positions to the literal vertices. Say variable A,; 
occurs unnegated in clauses l\ and I2 and negated in l 3 . Note that lj depend on 
the variable we are considering. For instance, in our example logical formula on 
p. [n| for variable A 1; li = 1, 1% = 3, and l 3 = 2. But for variable A 2 , h = 2, 
I2 = 3, and ^3 = 1. 

The preliminary literal vertex x\ is the intersection of the y = u(2l\ — 1) 
plane with the line connecting z l F and z l B . We do the same for the other positive 
occurrence of Xi and obtain the preliminary x\ . For the negative occurrence of 
Xi, we take the line connecting z l T and z B , intersect it with the y = u(2l^ — 1) 
plane, and obtain the preliminary x\. We join the preliminary x\ and bi 1 by a 
line d\ (this line lies in the y = u(2li — 1) plane). Do the analogue process for 
x\ and x\, obtaining d\ and d\. Later we will move the vertices x\, x\, x\ along 
their respective lines d\, d\, d\ a little out of polytope in order to turn them 
into extreme points. The lines dj will also be used for blocking conditions. 

Let H l be the plane that contains z % T and z F and the midpoint of the edge 
(z l L ,z B ) (Figure HJ). The only vertices above H l are x\, x\, x\, z\, and z B , 
and the only vertices on W are z % T and z F . This follows from the convexity of 
the current polytope. 

Let g\ (g l 2 ) be the line in the plane H l which is incident to z F and intersects 
the line d\ (d\). Note that this intersection point lies in the segment (x\,C2i 1 -i) 
(the line segment (a? 2 , C2i 2 -i)), thus in the interior of the constructed polytope. 
Analogously, let 173 be the line in the plane H l which is incident to z T and 
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y = it(2?i — 1) plane 





Figure 23: Construction of the literal vertices in the slice of the wedge 




Figure 24: Construction of H l and g\ 



intersects the line segment (x\, C2/3-1). It can be verified that if the roof height 
is small (z l L , z l R , z B ) is pierced by the g l j in its relative interior. This is another 
strict polynomial inequality in t roo f. It will be the planes H l and lines g l j (i = 
1, . . . , V) from which we make the visibility cones for the cupolas of variables 



Xi according to Theorem 2.11 



It is important to note right now that the non-blocking conditions are sat- 
isfied for this special position of the vertices. We do not want the tetrahedron 
{ztp, z % L , z Rl z l B ) and the triangle (x\, C2; 1 _2, c 2h) to intersect . From this we get 
strict polynomial inequalities on i roo f- They are satisfied for t roo f = since the 
y coordinates of the spine vertices q are smaller than 1/2. A suitable value of 
t roo f can be found solving the univariate inequality system we accumulated in 



our discussion (Lemma 1.4). It is easy to check that the sweeping conditions 
are also satisfied for the preliminary position of the points x\ , x\ , x\ . So far we 
have met two of the five required conditions to have a logical polytope. 
Stage 3: Literal vertices Now we put the final x l j (j = 1,2,3) a little outward 
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on line dj (Figure ^). A little for x\ and x\ means that the positive literal 
vertices lie in a plane parallel to the face (z R , z B , z\, z F ) very close to it. We 
treat x\ similarly. If the three literal vertices are moved a sufficiently small 
distance f literal the face lattice of what we get after taking the convex hull is as 



Figure 15 in all roofs. 

By construction Hi contains z F and z l Tl and the y = u(2j — 1) planes contain 
all literal vertices corresponding to clause j. This will become important for the 
visibility conditions (see Stage 5). Also, for small literal the non-blocking and 
sweeping conditions are satisfied. 

Although we do not have the blocking condition yet auxiliary lines can be 
set up: As above, let l\, I2, h be the clauses to which the literal vertex x\, 
x 2 , x\ belong. We made sure that the line segments (c^-i, x\) and (z F ,z B ) 
intersect in their respective relative interiors. Hence, by the construction of 
line g\, it is also pierced by (x\,C2i 1 -i) between zf and the face {z\,z % R ,z B ). 
(Analogously, (c2i % -\,x\) and (z F ,g 2 f~l (z l L , z R , z l B )) as well as (c2i 3 -\,x\) and 
(z T ,gl n (z l L , z R , z B )) intersect in their relative interiors). Later on this inter- 



section will evolve into the real blocking conditions using Theorem 2.11. 
Stage 4'- The perturbing the vertices on the spine of the wedge. We now perturb 
the points Cj on the spine of the wedge. Every even-indexed C21 is changed to lie 
on a parabola, and for the moment the odd-indexed vertices C2i-\ arc changed 
to lie on the line connecting 021-2 and C2i- The y coordinates of all points stay 
the same: 

C21 = Q(y - I) 2 • icvon, y, (y - l) 2 ■ t cv 

Note that by the 1/2 in the x coordinate, the points are moved into the 
polytope. The changes (parameter i evGn ) must be small enough that the convex 



hull now has the desired appearance (Figure 25) and the non-blocking conditions 



and the sweeping conditions arc still satisfied. Once more we appeal to Lemma 



1.4 The polynomials inequalities are now on the variable t cvcn and the sweeping 
and non-blocking were satisfied at t OVC n = 0. The reader should note that 
while the constructed vertices in the roofs do not change coordinates, dependent 
construction elements like the lines dj (connecting Xj and ca^-i) and (lying 
in H l and intersecting dj) change when the spine vertices move. However, the 
parameter i e ven has to be small enough that the preliminary blocking conditions 
are still met: g* still pierce the facet (z l L , z R , z B ) in its relative interior, and 
and d l j intersect in the interior of the polytope. 

Now we move the odd points C2/-1 beyond the face Gi = (021-2, C21, z^): to 
this end we choose a point pi beyond Gi and move to C21 - 1 +< dd (jpi — C21 - 1 ) . Such 
a point pi is easily found by taking a normal to G/ through its barycenter and 
moving out wards while staying beyond the face (note that this involves again 



Lemma 1.4, see the definition of beyond). The parameter t ^d is chosen small 
enough: Convexity and the correctness of the face lattice are easily achieved. 
Also the sweeping conditions are valid for slight moves. Keeping i oc id small also 
guarantees the non-blocking conditions: the tetrahedron (x\, 02^-2, 02^-1, C2i ± ) 
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Figure 25: Perturbation of the vertices on the spine 



is only slightly bigger than just the triangle [x\, 02^-2, C2h) which did not in- 
tersect the tetrahedron (z % Tl z l L , z R , z B ) (x 2 and x l 3 ). 

For the blocking conditions, let Xi be the jth logical variable in clause I. 
Note that now the line d* intersects the triangle (c2i-2, C2/-1, C21) in its relative 
interior. The lines are updated as the lines d* move. Since i dd is small, g^ 
still pierces the facet (z l L , z R ,z B ) in its relative interior, and g* and d* intersect 
in the interior of the polytope. Note that d* is still in the y = u(2l — 1) plane 
because the y coordinates of the spine vertices were conserved. 
Stage 5: Attaching the cupolas. It remains to construct all the cupolas. Over 
the facets (z l L ,z R , z B ) (i = 1, . . . , V) we construct cupolas using the planes Hi 
and sets of lines {g\,g2, gl}, and over the facets (C21-2, C2i-i, C21) (i — 1, ■ • ■ , C) 
we construct the clause cupolas using the y — g(2l — 1) planes and the sets 



of lines {d*- | X^s jth occurence is in clause I}. We invoke Theorem 2.11 and 
get the final polytope. By this construction, it is convex, has the correct face 
lattice, and the visibility conditions are satisfied. 

The reader will recall that g!j and dj intersect in the interior of the polytope. 
Say again variable Xj occurs unnegated in clauses l\ and I2 and negated in 



I3. By Theorem 2.11 5* pierces the skylight of the cupola corresponding to 
variable Xi and d* pierces the skylight corresponding to its clause lj. Hence, the 
tetrahedron spanned by Zp and the variable X^s skylight and the tetrahedron 
spanned by x\ (x 2 ) and clause Zi's skylight (fa's skylight) intersect in their 
interiors. Analogously, the tetrahedron spanned by z l T and the variable Xj's 
skylight and the tetrahedron spanned by x\ and clause l^'s skylight intersect in 
their interiors. These are exactly the blocking conditions. 

All other conditions concerned only points we constructed before, so they 
are still satisfied. The final polytope is therefore a logical polytope. □ 



Proof of Theorem 1.1: The problem is clearly in NP: checking whether a 
collection of tetrahedra is indeed a triangulation of the polytope P needs only 
a polynomial number of calculations. Every pair of tetrahedra is checked for 
proper intersection (in a common face or not at all), and the sum of the volumes 
equals the volume of P (computable for instance by the Delaunay triangulation 
of the polytope) . Also the size of triangulations of a given polytope is bounded 
by a polynomial in n of degree two (this follows from the well-known upper 
bound theorem, for details sec pl[). 
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By Lemma 3.3, from a given logical formula on V logical variables and C 
clauses, we can construct a logical polytope P of encoding length polynomial 
in V and C. Hence, by Lemma 3.2 there is a polynomial transformation that 
establishes the polynomial equivalence of a solution for the SAT problem and 
small triangulations of P. This completes the proof. □ 



4 Related Problems and Conclusions 



In this last section we discuss Corollary [L2 stated in the introduction and 
relate our main theorem to previously known results in the literature of optimal 
triangulations about the size of triangulations. We conclude with some remarks. 



Proof of Corollary \LQ: For part (1), consider the logical polytopes P we 
have constructed. If we could find the size s m i n for their minimal triangulations 
in polynomial time, either s m i n < K , in which case we have indeed a triangula- 
tion smaller or equal to K, or K < s mi „ in which case we can be sure there is 
no triangulation for the logical polytope of size K or less. This proves it must 
be hard to find a minimal triangulation in dimension 3. To extend to any other 
dimension simply note that pyramids over the logical polytopes have triangula- 
tions determined essentially by the logical polytope. In fact part (2) follows also 
from the pyramid construction because the logical polytope is a face. □ 

It is known that the sizes of triangulations for a d-dimensional polytope 
with n vertices lie between n — d and fd(dC(n + l.d + 1)) — d — 1, where 
fd(dC(n+ 1, d+1)) is the number of d-dimensional facets of a d+ 1 dimensional 
cyclic polytope with n+1 vertices [^lj . In particular, for 3-polytopes the possible 
number of tetrahedra ranges from n — 3 to — 2?i + 3. Both bounds arc known 
to be tight for three dimensions . It is also known that the size of a minimal 
triangulation of a convex 3-polytope must lie between n — 3 and 2n — 10, when 
n > 12 fl3|j . That the upper bound is tight was proved in |2q| using hyperbolic 
geometry. It is worth noticing at this point that the size of the constant K wc 
constructed in the previous section satisfies in fact n — 3 < K < In. Now wc 
discuss an interesting reason why the lower bound is strict: 

Wc say that a convex polytope is stacked if it has a triangulation whose dual 
graph is a tree (the dual graph of a simplicial complex is the graph that has one 
vertex for each maximal-dimensional simplex and two vertices are connected 
precisely when the corresponding simplices are adjacent via a common facet). 
The reader should be aware that in the literature the terminology stacked poly- 
tope is often restricted to simplicial polytopes. Here of course, we use it allowing 
that the stacking of simplices may give coplanarities. For example, any 3-cube 
or triangular prism is a stacked polytope under our definition. 

It turns out that a convex d-polytope P with n vertices has a triangulation 
of size n — d precisely when P is a stacked polytope (see |2l]]). It is natural 
to ask which polytopes are stacked. If it were NP-hard to recognize stacked 
polytopes then this would provide another proof that the problem of finding 
minimal triangulations is also in the same class. However we can prove: 
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Theorem 4.1 1. For a convex d-polytope P with n vertices there is a 
0(n d+2 ) algorithm to decided whether P is stacked (i.e., P has a tri- 
angulation with n — d maximal simplices). The algorithm uses only the 1- 
skeleton of the polytope. In particular, the size of a minimal triangulation 
of a stacked polytope does not depend on the particular coordinatization, 
but only depends on its face lattice. 

2. A convex ^-dimensional polytope P is stacked if and only if its graph does 
not contain as a minor the graph of an octahedron or a pentagonal prism. 

Proof: We need some definitions that have been introduced earlier in the graph 
theory literature Jl|, ||, We say a graph G is k-decomposable if G has k + 1 
or fewer vertices or there is a subset of vertices S of G with at most k vertices 
such that 1) S is a cut, i.e. G — S is disconnected, and 2) each of the connected 
components of G — S has the property that when the vertices of S are added 
back together with the complete graph on those vertices, the resulting graph is 
again fc-dccomposable. 

It was shown in Theorem 2.7 of [Q that the class of fc-decomposable graphs 
is the same as the class of partial fc-trees: a graph is a k-tree if it can be reduced 
to the complete graph Kk , by a finite sequence of removals of degree k vertices 
with completely connected neighbors (i.e. neighbors of the vertex induce a 
complete graph Kk). A partial k-tree is simply an edge-subgraph of a fc-tree. 

Now we claim that a d-dimensional convex polytope P is stacked, if and only 
if its 1-skeleton is a partial d-tree. Here is the proof: the "only if" implication is 
clear from the definition of stacked polytope. We can prove the "if" implication 
by induction on n. The theorem is trivial if n = d + 1 because then P is 
a simplex and its graph is a complete graph, thus is a d-tree. Assume then 
n > d + 1 and that the result is true for polytopes with fewer than n vertices. 
Remember that if 1-skeleton G(P) is a partial d-tree then it is d-decomposable. 
Thus there is a cut S of cardinality at most d. The set S must have in fact 
cardinality d because G(P) is d-connected by Balinski's theorem. If one has a 
vertex-cutset S of cardinality d in G(P), then the hypcrplane H(S) spanned by 
S intersects G(P) only in the vertices of G(P) H S and in no edges (otherwise 
S is not a cut). In conclusion, H(S) l~l P is a (d — l)-simplex and because G(P) 
is ^-decomposable we can apply induction hypothesis to prove the polytopes 
H(S) + (~l P and H(S)~ (1 P are d-decomposable, and thus they are partial d- 
trees with fewer vertices than P, so both polytopes are stacked. Finally, note 
that their stacked triangulations match well at the common boundary simplicial 
facet H (S) (~l P. proving that P is stacked. 

Arnborg et al. discuss in |2j an algorithm that, for fixed values of fc, decides 
whether a given graph is a partial fc-tree. The total running time is 0(n k+2 ). 
We also know that computing the 1-skeleton of a fc-polytope can be done with 
the same complexity. This, together with the above claim completes the proof 
of the first part of the theorem. 

For the second part we observe that partial fc-trees form a minor closed 
family. This means that the set of partial fc-trees is closed under taking edge- 
deletion or edge-contraction operations. The famous results of p4| imply that 
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they can be characterized by a finite set of forbidden minors. El-mallah and 
Colbourn jl2| proved that a graph is a planar 3-tree if and only if it has no 
minor isomorphic to the graph of an octahedron or a pentagonal prism. This 
fact together with our claim complete the proof of the second part. □ 

Remarks: 

1) The "coning" triangulation proposed in [ p"3"[ provides an algorithm which 
is polynomial on the number of vertices and gives a 2-approximation of the 
minimal triangulation as it produces a triangulation of size less than or equal 
to 2n - 7. 

2) Given a 3-dimensional convex polytope P, a proper subset S of tetrahcdra 
with vertices in vertices(P), and a positive integer K. Deciding whether there 
is a triangulation of P that uses K simpliccs from S can also be proved to be 
an NP-hard problem using the constructions we explained. 

3) It is interesting to note that the constructions presented in || prove 
also that covering a convex 3-polytope with tetrahcdra can be done with fewer 
pieces than triangulating. A cover is a collection of simplices whose union is the 
whole polytope, but the elements can intersect in their interiors. It is would be 
interesting to know what is the computational complexity of finding minimal 
simplicial covers. Another interesting question would be what is the complexity 
of deciding whether a triangulation of the boundary of convex non-simplicial 3- 
polytopc extends to a triangulation of the whole polytope without adding new 
interior points? If hard, this could be used to establish another proof our results. 
The curious reader can easily prove that already for a triangular prism not all 
triangulations of the boundary extend to a triangulation of the whole polytope. 
Even more interesting. The triangulations of the boundary of a 3-cube extend 
or not depending on the coordinates of its vertices. 
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